GaN:Gd: A superdilute ferromagnetic semiconductor with a Curie temperature above 300 K 
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We investigate the magnetic and magneto-optic properties of epitaxial GaN:Gd layers as a function of the ex- 
ternal magnetic field and temperature. An unprecedented magnetic moment is observed in this diluted magnetic 
semiconductor. The average value of the moment per Gd atom is found to be as high as 4000 fls as compared 
to its atomic moment of 8 (is- The long-range spin-polarization of the GaN matrix by Gd is also reflected in the 
circular polarization of magneto-photoluminescence measurements. Moreover, the materials system is found 
to be ferromagnetic above room temperature in the entire concentration range under investigation (7x lO'^ to 
2x lO'^ cm^^). We propose a phenomenological model to understand the macroscopic magnetic behavior of the 
system. Our study reveals a close connection between the observed ferromagnetism and the colossal magnetic 
moment of Gd. 

PACS numbers: 75.50.Pp, 76.30.Kg, 75.50.Dd,71.35.Ji 



Efficient spin injection into semiconductor structures is one 
' of the most important requirements to realize future spintronic 
. devices. A ferromagnetic material compatible with semicon- 

■ ductor growth which thus could easily be integrated into semi- 
conductor devices would be the best choice. Since the discov- 

' ery of a Curie temperature Tc as high as 110 K in (Ga,Mn)As 
. yl, III-V dilute magnetic semiconductors (DMS) have gained 
' a lot of attention. GaN-based DMS have the potential for ex- 
hibiting high temperature ferromagnetism which is essential 
for developing future spintronic devices |2]. Several transition 
metal (TM) doped GaN DMS have been investigated. Room- 

■ temperature ferromagnetism in these materials, however, is of- 
ten found to result from the formation of precipitates rather 
than from a homogeneous alloy |3]. The rare-earth (RE) el- 

, ements could be an interesting alternative to transition met- 
' als. RE atoms have partially filled / orbitals which carry mag- 
netic moments and may take part in magnetic coupling like 
in the case of TM with partially filled d orbitals. While it is 
expected that the magnetic coupling strength of / orbitals is 
much weaker than that of d orbitals due to the stronger local- 
ization of the / electrons, Gd is the only RE element which 
has both partially filled 4/ and 5d orbitals. Together these or- 
bitals can take part in a new coupling mechanism proceeding 

■ via intra-ion Af-5d exchange interaction followed by inter-ion 
, 5d-5d coupling mediated by charge carriers 01 • 

Here, we perform a systematic study of the magnetic prop- 
erties of the DMS GaN:Gd. The magnetization measure- 
ments reveal an unprecedented magnetic moment of up to 
4000 /Ib per Gd atom. This colossal magnetic moment can 
I be explained in terms of a long-range spin-polarization of the 
' GaN matrix by Gd. This conclusion is independently con- 
firmed by the pronounced changes of the circular polarization 
in magneto-photoluminescence measurements upon light Gd- 
doping. Moreover, all samples are found to be ferromagnetic 
above room temperature. A phenomenological model is de- 
veloped to explain the macroscopic magnetic behavior of the 
system, suggesting a close connection between the observed 
ferromagnetism and the colossal magnetic moment of Gd. 

The 400-700 nm thick GaN layers with a Gd concentra- 
tion ranging from 7x10^^ to 2x10'^ cm^^ were grown di- 
rectly on 6H-SiC(0001) substrates using reactive molecular- 



beam epitaxy (RMBE). The RMBE system (base pressure: 
2x10^^" Torr) is equipped with conventional effusion cells 
for Ga (purity: 7N) and Gd (purity: 4N) |5] and an unheated 
NH3 (purity: 7N) gas injector. A substrate temperature of 
810°C (our standard growth temperature for GaN) was used. 
The Gd/Ga flux ratio was changed in order to adjust the Gd 
concentration in the layers. Nucleation and growth was mon- 
itored in situ by reflection high-energy electron diffraction 
(RHEED). A spotty (1x1) pattern, reflecting a purely three- 
dimensional growth mode, was observed during nucleation of 
the layers. The pattern quickly became streaky, reflecting two- 
dimensional growth. 

The Gd concentration of the layers was determined by sec- 
ondary ion mass spectrometry using a CAMECA IMS 4f sys- 
tem, employing O2 primary ions with an impact energy of 
7 keV. Mass interference due to parasitic molecular ions was 
avoided by selecting precisely the masses of two Gd isotopes 
(155.922 and 157.924). The resulting ion rates were identical 
to within 1%. Measurements at different locations, each of 
which involved an area of 50x50 /xm^, deviated by no more 
than 5%. The concentration was calculated according to mea- 
surements of an ion-implanted standard (sample I), which was 
implanted with a Gd dose of lO'^ cm^^. The detection limit 
was found to be 2x lO'^ cm^^. 

Magnetization measurements up to 360 K were done in a 
Quantum Design superconducting quantum interference de- 
vice (SQUID) magnetometer with the sample held in a plas- 
tic straw provided by Quantum Design. The response of the 
magnetometer was calibrated with epitaxial Fe and MnAs lay- 
ers on GaAs for which the magnetization is accurately known. 
The magnetization of undoped GaN layers was found to be 
indistinguishable from that of bare SiC substrates. Magnetiza- 
tion loops were recorded at various temperatures for magnetic 
fields between ±50 kOe. The magnetic field was applied par- 
allel to the sample surface, i. e., perpendicular to the c-axis, 
for all measurements. All data presented here were coiTected 
for the diamagnetic background of the substrate. 

For magneto-PL, the sample was excited by a He-Cd ion 
laser at a wavelength of 325 nm. The laser power densities 
focused on the sample were in the range between 200 and 
250 Wcm^^. The experiments were caiTied out in the tem- 



2 




Flux ratio Jq^Jq^ 

FIG. 1: Gd concentration as measured by SIMS as a function of 
Gd/Ga flux ratio (solid squares). The solid line is a linear fit to the 
data representing sample C, E and F. The Gd concentration for sam- 
ple A, B and D (open squares) is obtained from the corresponding 
Gd/Ga flux ratio by extrapolation. The inset of the figure shows the 
SIMS profile obtained for sample C. 

perature range 5-100 K, in a continuous He-flow cryostat, and 
in magnetic fields up to 12 T using the backscattering Fara- 
day geometry. In all measurements, the magnetic field was 
oriented parallel to the c-axis of the samples. The magnetic 
field induced circular polarization of the (D^jX) transition was 
analyzed by passing the emitted light through a photoelastic 
modulator (PEM), a quarter-wave-plate, and a linear polarizer 

GaN:Gd layers A-G were grown with the Gd/Ga flux ratio 
increasing from A to G in alphabetic order. The concentration 
of Gd (A^Gd) in samples C, E, F and G is measured by sec- 
ondary ion mass spectrometry (SIMS). In Fig.^ A^cd as mea- 
sured by SIMS is plotted as a function of the Gd/Ga flux ratio 
([>. Samples C, E and F (solid squares) are lying on a straight 
line with a slope of unity indicating a linear dependence of 
Gd incorporation on (j). The Gd concentration of sample G is 
smaller than expected from this linear dependence. In fact, we 
have observed a strongly faceted surface for sample G, indicat- 
ing that Gd is affecting the growth mode at these high concen- 
trations. Ncd for sample A, B and D (open squares) is obtained 
by Unearly extrapolating the curve passing through samples C, 
E and F. The inset of the figure shows the SIMS profile ob- 
tained for sample C. The concentration of Gd is found to be 
constant over the entire depth for sample C and all other sam- 
ples investigated by SIMS, ruling out any accumulation of Gd 
on the surface during growth. All samples were subject to an 
extensive investigation by high-resolution x-ray diffraction in 
a wide angular range. No reflections related to a secondary 
phase was detected. Furthermore, we have investigated one 
sample (C) by transmission electron microscopy. No clusters 
or precipitates were observed. Electrically, all Gd doped sam- 
ples are found to be highly resistive. 

Figure |2] (a) shows the magnetization loops obtained for 
sample C at 2 K and 300 K. At both temperatures, the magneti- 
zation saturates at high magnetic fields (unlike superparamag- 
netic materials) and exhibits a clear hysteresis at lower fields. 
These two features clearly indicate a ferromagnetic behavior 
We have observed qualitatively similar hysteresis loops for all 



samples. Figure |2](b) and (c) show the temperature depen- 
dence of the magnetization under a magnetic field of 100 Oe 
for samples A and C, respectively. Prior to measuring the tem- 
perature dependence of the magnetization, the sample is first 
cooled from room temperature to 2 K either under a saturation 
field of 20 kOe (Field cooled: EC) or at zero field (zero field 
cooled: ZFC). It is clear from these figures that the FC and 
ZFC curves are qualitatively similar for both samples featuring 
a double step like structure below 100 K in the FC curves. In 
the case of sample C, the two curves remain separated through- 
out the entire temperature range (2-360 K), while they coin- 
cide at around 360 K for sample A. The separation between 
the FC and ZFC curves indicates a hysteretic behavior which 
is consistent with our previous observation from Fig. |2l (a). 
The two curves coincide at the Curie temperature T^, when the 
hysteresis disappears. Clearly, 7^ is around 360 K for sample 
A while it is much larger for sample C. The separation be- 
tween the FC and ZFC curves measured at 360 K is plotted as 
a function of Gd concentration in the inset of Fig.|2l(c), and is 
seen to increase with increasing Gd concentration revealing a 
shift of Tc towards higher temperatures. We have also inves- 
tigated a Gd-implanted GaN sample (sample I) with a higher 
Gd concentration than we could reach by incorporation during 
growth. This sample also exhibits ferromagnetism at and well 
above room- temperature. These observations are consistent 
with the results obtained by Teraguchi et al, who observed a 
Curie temperature larger than 400 K in Gao 94Gdo oeN alloys 
16]. 

The effective magnetic moment per Gd atom p^f^ can be ob- 
tained from the value of the saturation magnetization {pes 
= Ms/Nq^). In sample C, peff is found to be 935 and 737 jj.^ at 
2 K and 300 K, respectively. These values are about two order 
of magnitudes larger than the pure moment of Gd. In Fig.|3l 
Peff obtained at different temperatures is plotted as a function 
of A^Gd- Data obtained from sample I are also included in this 
figure (solid squares). Both at 2 K [Fig.|3l(a)] and 300 K [Fig. 
|3l(b)], psff is found to be extraordinarily large, particularly for 
low Gd concentrations. At high Gd concentrations, peff is seen 
to saturate. It is interesting to note that at 300 K, the satura- 
tion value is close to that of the atomic moment of Gd (8 JJ-b), 
while at 2 K, it is clearly higher. This is true even for sam- 
ple I. Such a colossal moment can be explained in terms of a 
very effective spin-polarization of the GaN matrix by the Gd 
atoms. In fact, in certain diluted metallic alloys, the solute 
atoms exhibit an average magnetic moment larger than their 
atomic value. This effect is called giant magnetic moment. 
There are several reports on a giant magnetic moment of Fe, 
Mn and Co when they are either dissolved in or residing on 
the surface of Pd j^laBl- This effect is explained in terms 
of a spin polarization of the surrounding Pd atoms. Very re- 
cently, a giant magnetic moment was observed in Co doped 
Sn02_5 IIOII . a DMS. High-temperature ferromagnetism and 
an unusually large magnetic moment were observed also in 
Fe-doped Sn02, which was argued to originate from ferrimag- 
netic coupling mediated by electrons trapped on bridging F- 
centers 1 1 y . 

The above conclusion of a long-range spin-polarization of 
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FIG. 2: (a) Magnetization loops obtained for sample C at 2 K (solid 
line) and 300 K (dotted line). The inset of (a) shows the loops at small 
fields. Temperature dependence of magnetization at field-cooled (FC) 
and zero-field-cooled (ZFC) conditions at a magnetic field of 100 Oe 
for (b) sample A and (c) sample C. The difference between FC and 
ZFC magnetization measured at 360 K is plotted as a function of A'c^j 
in the inset of (c). The solid curve through the data is guide to the 
eye. 



the GaN matrix by the Gd atoms is independently confirmed 
by magneto-PL. In the absence of an external magnetic field, 
the PL spectra of all samples are characteristic for undoped 
epitaxial GaN layers in that they are dominated by a neutral- 
donor bound exciton (D'',X) transition at ■~3.458 eV. The 
lower energy of this transition when compared to homoepitax- 
ial GaN is consistent with the tensile in-plane strain in these 
layers of 0.15% fli"). In all cases, the donor responsible for 
this transition is O with a concentration of about lO'*^ cm^^ as 
measured by SIMS. However, clear differences between un- 
doped and Gd-doped GaN layers are evident when analyzing 
the circular polarization of the PL in the presence of a mag- 
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FIG. 3: Magnetic moment per Gd atom (Peif) as a function of Gd 
concentration at (a) 2 K and (b) 300 K. The saturation magnetization 
(Ms) as a function of Gd concentration at 2 and 300 K is also shown in 
the respective insets. Open squares are representing the experimental 
data for sample A-G and solid squares for sample I. Solid lines are 
the theoretical fit obtained from our model (discussed in the text). 



netic field. Figure@]shows the PL spectra for a GaN reference 
sample [Fig. 0](a)] and sample C [Fig. |4](b)] under a mag- 
netic field of B = 10 T. The transition is polarized for both 
samples, which is evident from the difference in intensities of 
the two circularly polarized (7+ (open squares) and a (solid 
squares) components. Most importantly, the polarization has 
opposite sign for sample C when compared to the reference 
sample, which clearly indicates that the excitonic ground state 
in sample C is subjected to an internal magnetic field which is 
effectively stronger than the external field and is acting in the 
opposite direction. Considering that the mean Gd-Gd separa- 
tion in sample C is 25 nm, and that the average distance of a 
(D'',X) site from a Gd atom will thus be as large as 12 nm iTsll . 
this finding directly confirms our previous conclusion that the 
Gd atoms are introducing a long range spin polarization in the 
GaN matrix. We also observe that the relative change of po- 
larization with respect to the reference sample increases with 
the Gd concentration, as intuitively expected. These results, 
together with a quantitative analysis of the dependence of the 
polarization on the magnetic field, will be reported in detail 
elsewhere 1 14]. 

In order to get a quantitative understanding of the range of 
the spin polarization, we have developed a phenomenological 
model as explained in the following. The polarization of the 
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FIG. 4: Circularly polarized PL spectra for (a) the reference GaN 
sample and (b) sample C, measured at 6 K under a magnetic field of 
6 = 10 T in Faraday configuration (B \\ c-axis). 

GaN matrix by the randomly positioned Gd atoms is described 
as a rigid sphere of influence around each Gd atom, meaning 
that all the matrix atoms within the sphere are polarized by 
an equal amount whereas matrix atoms falling outside of this 
sphere are not affected. Let us associate an induced moment of 
po with each of the matrix atoms lying in the region occupied 
by one sphere. As one should expect an increase of the polar- 
ization of matrix atoms if they belong to a region where two (or 
more) spheres of influence overlap, we attribute an additional 
moment npito matrix atoms occupying sites in regions where 
n (n ~ 2,. ..Nad) spheres overlap. Since their radius r is pre- 
sumably much larger than the lattice spacing in GaN, we fur- 
ther assume that the spheres are randomly arranged in a three- 
dimensional continuum (continuum percolation). Within this 
framework, the saturation magnetization can be expressed as 

Ms = f GdAfed + PovNo + piNq nv„ (1) 

n=2 

where A'o is the concentration of matrix atoms per unit volume, 
V is the volume of each sphere, 1 — exp(— vA^od) is the 
volume fraction occupied by the spheres, and 



is the volume fraction of the regions contained within n 
spheres. The average magnetic moment per Gd atom /^gg- is 
then obtained as 



+ h-(/'o + /^iA^Gdv)e-"'^«] (3) 

-'VGd 

At low Gd concentrations, most of the spheres are well sep- 
arated and peff has its maximum value. As the concentration 
of Gd is increased, more and more spheres overlap and /?eff de- 
creases. Finally, at a very high Gd concentration the entire ma- 
trix becomes polarized, and /7eff approaches saturation. Note 
that the value of saturation is larger than the magnetic moment 
of bare Gd atoms by an amount of p\Nqv. We use Eqn.|3lto fit 
our experimental data with pGd = 8 /^b and po, pi and r as free 
parameters. The agreement is quite satisfactory as shown in 
Fig.|3] The fit yields po=l.l x IO^Vb,/?! = 1-0 x IO^Vb, 
and r = 33 nm at 2 K and po = 8.4 x lO"'* ^b, pi ~ 0, and 
r = 28 nm at 300 K. The finite value of pi at 2 K explains why 
Peff saturates at a value which is still higher than the atomic 
moment of Gd. Most important, however, is that the value for 
r is sufficiently large to account for the strong effect of Gd on 
the exci tonic ground state of the GaN matrix (cf. Fig.E}. 

The second remarkable property of our GaN:Gd layers is 
the high-temperature ferromagnetism observed in all samples 
(well above room temperature, compared to a Tc of 289 K for 
bulk Gd). We believe that the ferromagnetism and the colossal 
moment of Gd observed in these samples are in fact closely 
related. An overlap of the spheres of influence establishes a 
(long-range) coupling between the individual spheres. Within 
the framework of percolation theory, ferromagnetism is ex- 
pected to occur at the percolation threshold, when an "infi- 
nite cluster" spanning macroscopic regions of the sample is 
formed. The percolation threshold is reached for our model 
at V = 0.28955 [15, J_6]. With increasing Gd concentration, 
we would thus expect a phase transition from paramagnetic to 
ferromagnetic behavior. Furthermore, Tc which depends upon 
the strength of the overlap is expected to increase with the Gd 
concentration. This is clearly consistent with the results shown 
in the inset of Fig.|2](c). 

Unfortunately, the well-established percolation formalism 
cannot be straightforwardly adopted to quantitatively explain 
the current problem. It is intuitively clear that the spheres of 
influence are not hard but soft, i. e., the polarization of the 
matrix induced by the Gd atoms must decay with increasing 
distance. The onset of ferromagnetic order now depends upon 
the precise shape of this polarization cloud as well as on the 
overlap of two or more of these clouds. A quantitative pre- 
diction of the onset of ferromagnetism in this situation thus 
requires a detailed understanding of the nature of the ferro- 
magnetic coupling in this material. 

Regarding the microscopic origin of this coupling, it is clear 
that it cannot be explained simply in terms of direct, double 
or superexchange between Gd atoms since the average Gd- 
Gd distance is too large for such a coupling to exist. Further- 
more, all samples are found to be electrically highly resistive, 
ruling out free-carrier mediated RKKY (Rudermann-Kittel- 
Kasuya-Yosida) type long-range coupling. Finally, also very 
recent proposals aimed at explaining a giant magnetic moment 
and/or high-temperature ferromagnetism in various material 
systems llZl la. ISl Uil ll 111 do not apply to the cuiTent require- 
ment of a long-range coupling. An actual understanding of 
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the phenomenon observed in the present paper will require de- 
tailed ab initio studies, particularly of the f-d coupling be- 
tween Gd and Ga. Considering the long spatial range of the 
coupling evident from our experiments, such calculations are 
computationally very demanding and are therefore beyond the 
scope of the present paper An important result for any theoret- 
ical approach was recently obtained by Hashimoto et ai, who 
showed that Gd occupies primarily Ga sites in GaN even at a 
concentration as high as 6% (two orders of magnitude higher 
than the highest concentration investigated here) JitIi . 

Our study has revealed an extraordinarily large magnetic 
moment of Gd in GaN which can be explained in terms of 
a long-range spin-polarization of the GaN matrix by the Gd 
atoms. The material is found to exhibit ferromagnetism above 
room temperature even with a Gd concentration less than 
lO'^ cm^^. This finding offers an exciting opportunity, since 
GaN:Gd may be easily doped with donors (acceptors) with a 
concentration exceeding that of Gd to generate spin polarized 
electrons (holes) in the conduction band (valence band). Gd- 
doped GaN with its 7^ above room temperature might thus be a 
very attractive candidate as a source of spin-polarized carriers 
for future semiconductor-based spintronics [18']. 
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